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By properly driving a qubit-resonator system in the strong dispersive regime, we implement an 
“impedance-matched” A system in the dressed states, where a resonant single photon deterministi¬ 
cally induces a Raman transition and excites the qubit. Combining this effect and a fast dispersive 
readout of the qubit, we realize a detector of itinerant microwave photons. We theoretically analyze 
the single-photon response of the A system and evaluate its performance as a detector. We achieve 
a high detection efficiency close to unity without relying on precise temporal control of the input 
pulse shape and under a conservative estimate of the system parameters. The detector can also be 
reset quickly by applying microwave pulses, which allows a short dead time and a high repetition 
rate. 

PACS numbers: 42.50.Pq 03.67.Lx 85.25.Cp 

I. INTRODUCTION 

Extensive efforts have been made in a varieW of physical systems to realize strong coupling between a single quantum 
emitter and a one-dimensional photon field |ll-ll0l|. In such one-dimensional optical systems, interaction between an 
emitter and a photon is enhanced drastically due to the destructive interference between the incident field and the 
radiation from the emitter. This opens the possibility for deterministic control of quantum systems by individual 
photons. 

In particular, in a A-type three-level system that has identical radiative decay rates from the top level and is coupled 
to a semi-infinite one-dimensional field, a resonant single photon deterministically induces a Raman transition in the 
A system and switches its electronic state PH. Recently, we implemented such a A system by utilizing the 
dressed states of a driven circuit-quantum-electrodynamics (circuit-QED) system. Applying continuous microwave to 
the A system, we confirmed that the field amplitude vanishes completely upon its reflection (perfect absorption by 
“impedance matching”), and that input photons are down-converted by the A system PP. This indicates highly 
efficient switching of the A system induced by individual microwave photons. Deterministic switching of a A system 
has also been demonstrated recently in the visible-light domain by using a spherical cavity and an atom P, P] . 

Exploiting the large transition dipole of superconducting qubits, circuit QED enables various microwave quantum- 
optical phenomena that have not been reached by quantum optics in the visible domain; The examples include direct 
and strong coupling of qubits to itinerant waveguide photons |20l422l| and various effects in the strong-dispersive regime 
of the qubit-cavity interaction [23l - l^ . However, an apparent shortcoming of microwave quantum optics has been the 
lack of efficient single photon detectors: The main difficulty is in the energy scale of single microwave photons orders 
of magnitude smaller than that of visible or infrared photons. 

In recent pro posals and an experiment, a current-biased Josephson junction was used as a microwave photon 
detector [261 - 1^ : a metastable two-level system formed in the tilted washboard potential resonantly absorbs an 
incoming photon and switches into the voltage state of the junction to make a “click”. While the successful operation 
was demonstrated P|, the dark count and the dead time after the detection are remaining issues. In this study, 
we theoretically analyze the single-photon response of an impedance-matched A system realized in a driven qubit- 
resonator system and demonstrates its excellent performance as a single photon detector in the microwave domain. 
The detector has a low dark count rate and attains a high detection efficiency close to unity. Moreover, the system 
is switched to a discrete and definite quantum state upon detection, which enables us to reset the system quickly 
and therefore to shorten the dead time of the detector. In comparison with the works that demonstrate deterministic 
capture of propagating microwave pulses in a harmonic oscillator mode (29l - [3^ , our scheme based on a A system has 
a merit that it is free from precise temporal control of the input pulse shape and/or the system parameters: high 
switching efficiency close to unity can be attained as long as the input photon has a narrower bandwidth than the 
linewidth of the A system. 

The rest of this paper is organized as follows. In Sec. [Hi we theoretically describe the driven qubit-resonator system 
and derive basic equations for the analysis. In Sec. uni we explain two distinct functionalities of the device. In the 
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FIG. 1: Schematic of the setup. A superconducting qubit (two-level atom) is coupled dispersively to a resonator. The 
resonator is coupled to a semi-infinite waveguide (WGl), through which a signal photon is input. The qubit is coupled to 
another waveguide (WG2), through which a drive pulse is applied. 


“A” mode, which is realized under a proper qubit drive, a single incident photon excites the qubit deterministically. 
In the “I” mode, which is realized when the qubit drive is off, the dispersive readout of the qubit state is possible. 
In Sec. lYl we present our single-photon detection scheme, which is composed of three stages: capture, readout, and 
reset. In Sec. El we discuss the capture stage and show that a high detection efficiency close to unity is achievable 
with practical parameter values. In Sec. ED we discuss the reset stage and show that the present system can be 
initialized within a few hundred nanoseconds. Section IVIII summarizes this study. 


II. SYSTEM 

Here, we theoretically describe the device considered in this study [H, H^. Its schematic is shown in Fig. [1] A 
superconducting qubit, which can be regarded as a two-level atom, is dispersively coupled to a transmission line 
resonator. The resonator is further coupled to a semi-infinite waveguide (WGl), through which a signal photon pulse 
to be detected is input. Our objective is to determine whether the signal pulse contains a photon or not. We also use 
this port for readout of the qubit and for resetting the system. Through another waveguide (WG2), we apply drive 
pulses to the qubit in order to engineer the dressed states of the qubit-resonator system. 


A. Hamiltonian 


Setting h = V = where v is the microwave velocity in the waveguides, the Hamiltonian of the overall system is 
written as 


TL = T-Lsys + 'Hda 


'Hsys = uJra^aaa^ -I- [ujg + {wr — 2x)a)a\ crV, 
■Hdamp = dk ka\ak + \/ k'/ 2'n:{a^ak + a\a) 


dk 


kh\hk + yjx' l2'K{cF%k + h\c 


( 1 ) 

( 2 ) 

( 3 ) 


where a (a) is the annihilation operator for the resonator (qubit) and at (bk) is the annihilation operator for microwave 
photon propagating in WGl (WG2) with wave number k. iOr (j-Oq) is the resonance frequency of the resonator (qubit), 
X is the dispersive frequency shift, and k' ( 7 ') is the radiative decay rate of the resonator (qubit) into WGl (WG2). 
The parameter values are listed in Table ID 

Three comments are in order regarding this Hamiltonian, (i) The qubit-resonator system is described by the 
Jaynes-Cummings Hamiltonian, Hsys = -I- WgCrV -I- g{a^a- + a^a), where uJr and ujq are the bare frequencies of 

the resonator and qubit and g is their coupling. In the dispersive regime {\iUr — 0 Jq\ ^ g), Hsys is recast into a diagonal 
form of Eq. ©■ The renormalized frequencies are given hy cOr = tOr + X and LOq = ujq — x, where x = g'^li^r — Wg). 
(ii) Although omitted for simplicity, the resonator and qubit are subject to nonradiative decay to the environment. 
We denote the total decay rate of the resonator (qubit) by k ( 7 ). (iii) In practice, we switch to the rotating frame 
to remove the natural phase factors. We subtract Ho = uJsa^a + oJdO'^cr from Hsys of Eq. ([2]), where uJs (wd) is the 
central frequency of the signal photon (drive pulse). 
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OJq 

UJr 

X 

27r X 5 GHz 

27r X 10 GHz 

27r X 40 MHz 

K 

27r X 20 MHz 

k! 

27r X 20 MHz 

7 

27r X 0.1 MHz 

7' 

27r X 0.1 kHz 


ujq — 2'it X 70 MHz 


TABLE I: Parameter values for ujq (qubit frequency), Ur (resonator frequency), x (dispersive shift), k, (total decay rate of the 
resonator), k' (radiative decay rate of the resonator to WGl), 7 (total decay rate of the qubit), and 7 ' (radiative decay rate of 
the qubit to WG2), and ojd (qubit drive frequency). 


B. Heisenberg equations 


We introduce the real-space representation of the field operator for WGl by 


1 



dke'-’^'^ak. 


(4) 


In this representation, the r < 0 (r > 0) region corresponds to the incoming (outgoing) field. The field operator br 
for WG2 is defined similarly. From T^damp of Eq. ([3]), we can derive the following input-output relations: 

dr{t) = dr-t{0) — — r)9{r)6{t — r), (5) 

br(t) = br-tiO) - - r)9{r)0{t - r). (6) 

The conventional input and output field operators are defined at r = ±0 by ain(t) = d-o(t) and aout{t) = d^o{t). 
bin{t) and 6out(t) are defined similarly. From Eqs. ([I|), ([5]) and ([6|), the Heisenberg equation for any system operator 
S (composed of ct, a and their conjugates) is given by 

— S = ilT-Lsys, S] +—{2a^Sa — o'aS — Sa^a) + ^{2a^Sa — a^aS — Sa^a) 

+ i'/idal^{t)[a, S] + iVn'la^ S]ain{t) + i-J^b\^{t)[a,S] + S]bin{t). (7) 


C. Microwave response 


In our setup, a signal photon is input through WGl and a classical drive pulse is applied through WG2. We denote 
the wavefunction of the signal photon by fs{t) and the amplitude of the drive pulse by fdit). Note that fs{t) is 
normalized as f dt\fs{t)\‘^ = 1. Extension to the two-photon input through WGl is straightforward (see Appendix El). 

Analysis of the microwave response to a single photon can be simplified by replacing the single-photon state |1) 
with a coherent state |a), performing perturbation calculation with respect to a, and picking up the relevant terms 
afterwards [s^. Therefore, we investigate a situation in which two classical pulses, afs{t) and fd{t), are applied through 
WGl and WG2, respectively. Since a classical pulse (coherent state) is an eigenstate of an input field operator, we 
can rigorously replace ain(t) and &in(t) in Eq. (O with afs(t) and fd(t), respectively. Then, the expectation value of 
an operator S evolves as 

^(S')c = i{[Hsys, *S'])c + K{{a'‘Sa)c - {a^aS)c/2 - {Sa^a)c/2) + 7 ((cr'l'S'cr)c - {(J^aS)c/2 - (5'crV)c/2) 

-b f*{t){[a, 5'])c + iy/n'afs{t){[o', S])c + fd{t){V^, S])c+ iy/^S])c. ( 8 ) 

We expand {S)c as {S)c = Then, (S')™" evolves as 



= ^([iHsys, s])™" + K((a+Sa)™" - (ataS)™"/2 - (Sata)™"/2) + 7((atSa)™” - (aVS)™"/2 - (SaV)™"/2) 
+ *v^/;(t)([a, s])(™-i)"+^y;d^(^)([a^ s])r("-')+ s])r +5])r • ( 9 ) 


with the convention that (S)™" 


0 for negative m or n. 
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Our objective is to evaluate the expectation value for a single-photon input, (S')* 
is given by 


As discussed in Appendix it 


(S)* = (S)°° + (S)“. 


( 10 ) 


Therefore, we solve the simultaneous differential equations for (S))!°, (S))l^, (S)),° and (S)),^ to evaluate (S)*- 


III. I AND A MODES 

As discussed in Refs. [H, [l^, for a drive field with properly chosen power and frequency, the dressed states of a 
qubit-resonator system constitutes an “impedance-matched” A system, which absorbs an input photon with a high 
efficiency close to unity. To observe this, we first consider a case with a continuous drive in this section. The drive 
field fd{t) is given by 


m = ( 11 ) 

where is the drive frequency and fid is the drive amplitude expressed in terms of the Rabi frequency. In the frame 
rotating at ujd, the Hamiltonian for the driven qubit-resonator system is written as 

'Hays+dr = oJra^aaa"^ + [{ojq — ujd) + {ojr — 2x)a'''a] ctV -|- Od(cr'^ -I- a). (12) 

We label the state vectors of the system with |g, n), where q{= g, e) denotes the qubit state and n(= 0,1, • • •) denotes 
the resonator photon number. Throughout this study, only the lowest four states (|(?,0), |e,0), \g,l) and |e, 1)) are 
relevant. When the drive is off {fid = 0), their energies are respectively given by <jj|g,o) = Oj ‘^le.o) = ^q—^d, <-^\g,i) = 
and <jj|e,i) = ojq — bjd -\- idr — 2^. We choose the drive frequency Wd to satisfy ujq — 2x < uJd < ojq. Then, the level 
structure becomes nested, i.e., u;|g o) < <^|e,o> < Pig-ISKa)]. When the drive is on {fid > 0), |g, 0) and 

|e, 0) (|g, 1) and |e,l)) are mixed to form the dressed states |1) and |2) (|3) and |4)), where the dressed states are 
labeled from the lowest in energy [Fig. [2Db)]. The radiative decay rate from |i) to |j) emitting a photon into WGl is 
given by = K'|(i|af|j)p. 

Figure |2l[c) plots as functions of the drive amplitude. The drive frequency is fixed at the value in Table H] 
throughout this study. When the drive is off, 1532 = 154i = k! and Kai = 1542 = 0. Namely, the radiative decay occurs 
only in one direction as depicted in Fig. UDa). We refer to this as the “1” mode of the qubit-resonator system. In 
this mode, the microwave transition occurs conserving the qubit state, and the transition frequency differs by 2 x 
depending on the qubit state. We can use this mode for the dispersive readout of the qubit state. 

In contrast, at a proper drive power in Fig. |2l[c)], the four decay rates become identical, i.e., ksi = 1532 = 

1541 = 1542 = n'/2. Then, the three levels |1), |2) and |u) (m = 3 or 4) function as an impedance-matched A system: a 
single photon resonant to the |1) —^ |1I) transition deterministically induces the Raman transition of |1) —^ \u) —>■ |2) 
and switches the electronic state. We refer to this as the “A” mode of the qubit-resonator system [Fig. |2][b)]. The 
underlying physical mechanism for this phenomenon is the destructive interference between the input and elastically 
scattered photons [H- Figure[2Dd) shows the reflection coefficient |r| of a weak continuous signal field applied through 
WGl, as a function of the drive amplitude fid and the signal frequency w*. We observe two impedance-matching spots 
(|r| ~ 0) at {fld,u!s) — and (H™^,cti 4 i), where Wy is the transition frequency between dressed states |i) 

and \j). We have two comments regarding Fig.[2Id). (i) We observe that w* of the upper (lower) spot deviates slightly 
from W 41 (w 3 i). This is due to the small energy difference between levels |3) and |4) that is comparable to the linewidth 
of the A system (^ n'/2.). In this case, weak elastic scatteriiig induced by level |3) (|4)) is not negligible, (ii) In our 
former experiment, fid of the two spots differed considerably (iTI j. That was because we used a relatively strong signal 
field to improve the signal to noise ratio. Here, we consider the weak-field limit to discuss the single-photon response. 


IV. DETECTION SCHEME 

As we observed in the previous section, the present qubit-resonator system has two distinct functionalities. In the 
A mode, the system captures a single photon nearly deterministically and makes a transition to the qubit excited 
state. In the I mode, the system functions as an oscillator whose resonance frequency depends on the qubit state. We 
can read out the qubit state efficiently with a single probe pulse. These two modes can be switched adiabatically by 
a temporally smooth qubit drive pulse. 
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FIG. 2: (a) I mode of the qubit-resonator system. This is attained when the qubit drive is off. Arrows indicate the directions 

of radiative decays, (b) A mode of the qubit-resonator system. The four decay rates are identical. This is attained when a 
proper drive field is applied to the qubit, (c) Radiative decay rates K31, K32, K41 and K42 as functions of the drive amplitude. 
The four rates become identical at n™^/27r = 13.2 MHz. (d) Reflection coefficient |r| of a weak signal field applied through 
WGl. White dashed lines show a; 3 i and a; 4 i. Impedance matching (|r| ~ 0) occurs at {Q,d,u}p) ~ (f2™^,a53i) and (fI™^,0J4i). 


waveguide 1 




capture 


readout 


reset 


FIG. 3: Pulse sequence for single-photon detection, (i) Gapture stage. We simultaneously input a signal photon through WGl 
and a classical drive pulse through WG2. (ii) Readout stage. We apply a classical readout pulse through WGl and switch off 
a drive pulse through WG2. (iii) Reset stage. We simultaneously apply a classical reset pulse through WGl and a classical 
drive pulse through WG2. The carrier frequencies of the three pulses applied through WGl are different from each other. 


In the single-photon detection considered here, it is assumed that the arrival time and the length of the signal pulse 
are known in advance: we determine whether the pulse contains a single photon or not by this measurement. Our 
single-photon detection proceeds in the following three steps, (i) Initially, the qubit-resonator system is in its ground 
state |5,0) and a signal pulse is input through WGl. Through WG2, we apply a drive pulse that covers the signal 
pulse in time, and switch the system to the A mode during this period. At the end of this stage, the photon number 
in the signal pulse (0 or 1) is mapped onto the qubit state {g or e). We refer to this as the “capture” stage and discuss 
in detail in Sec. El (ii) Next, we perform dispersive readout of the qubit. In order to prevent the |2) —>• |m) —>• |1) 
transition (m = 3 or 4) induced by the probe, we turn off the qubit drive and keep the system in the I mode during 
this process. We apply a classical probe pulse to the resonator and measure the phase shift of the reflected pulse. We 
refer to this as the “readout” stage. We do not touch this stage in this paper, since dispersive readout of the qubit 
has been discussed in prior works: High-fidelity single-shot readout has been achieved with the Josephson parametric 
amplifier [s^ [s^ and the parametric phase-locked oscillator [s^. (iii) The qubit-resonator system is initialized (i.e., 
decays to its ground state \g, 0)) automatically by the natural qubit decay. However, in order to shorten the dead time 
of the detector, we artificially reset the system by applying microwave pulses. For this purpose, we use the inverse 
process of stage (i): We apply a drive pulse through WG2 to form a A system and a reset pulse through WGl. Note 
that the reset pulse can be a classical and strong one in contrast with the capture stage. We refer to this as the 
“reset” stage and discuss in detail in Sec. ED 


V. CAPTURE OF SIGNAL PHOTON 


In this section, we present the numerical results for the capture stage of Fig. [31 where a signal pulse to be detected 
and a qubit drive pulse are input simultaneously. We assume a Gaussian pulse with length I for the signal pulse: 


fs{t) 


81 n 2 y^'^ 

ttP ) 




(13) 


















6 


w/2 (3/ wl2 



FIG. 4: Envelopes of the signal photon and the qubit drive pulse. We use a smooth drive pulse in order to switch the I and 
A modes adiabatically. 


Note that fg is normalized as J dt\fs{t)\'^ = 1. On the other hand, we assume the following pulse shape for the qubit 
drive: 


fd{t) 


^ /1 i\t\ < aim 

VY |2-(I‘I-W2)V(-/2)^ (|i| > ;3//2) ■ 


(14) 


Namely, a square pulse with length fil is smoothened by Gaussian functions with width w. We smoothen the drive 
pulses in order to switch the I and A modes adiabatically and to suppress unwanted qubit excitations. A real positive 
constant /3 is chosen so that fd covers fs in time. For Gaussian signal pulses, we set /? = 2. The envelope functions of 
the signal photon and the drive pulse are drawn in Fig. |4l In the capture stage, we set the drive amplitude at 
indicated in Fig. He). 


A. Zero-photon input 

Here, we discuss the case in which the signal pulse contains no photons. Namely, we observe the dynamics induced 
solely by the drive pulse applied through WG2. In Fig.jSJ we plot the time evolution of the qubit excitation probability, 
Po{t) = m o). (The subscript 0 represents the photon number in the signal pulse.) For reference, we also plot the 
excitation probability Po(0 under the adiabatic approximation, which is obtained as follows. The system Hamiltonian 
including the qubit drive is written, in the rotating frame, as 'Hsys+dr = + [(w^ — w^) + (w^ — 2x)a^a] ( 7 V + 

•YY|/d(t)|(o'^ +cr). Its lowest two eigenstates |1) and |2) are given by |1) = cos0\g, 0) — sindje, 0) and |2) = sind|(/, 0) + 
cos0|e, 0), where 


9{t) = - arctan 


iOq- LOd ) 


(15) 


Under the adiabatic assumption, the system always stays in its ground state, |1). The excitation probability, Po(0 = 
|(e, 0|l)p, is therefore given by 


Po(0 = sin^^'W- (16) 

We observe in Fig.[5]that po(l) is at most a few percents and roughly agrees with The deviation between them 

originates in the nonadiabatic transition, which is suppressed by increasing w. Immediately after the drive is switched 
off, po(t) comes back close to zero and becomes nearly independent of time. Hereafter, we set tc = 30 ns and define the 
detection probability Pq as the qubit excitation probability at the end of this stage. Namely, Pq = poitf), where we 
set tf = pi/2 + 50 ns. Under this choice of w, the detection probability is negligibly small for the zero-photon input. 
This means that the dark counts are suppressed nearly completely in the present scheme. For example, Pn ~ 10“^ 
for I = 200 ns and /3 = 2. 


B. One-photon inpnt 

Here, we discuss the case in which the signal pulse contains one photon. We denote the qubit excitation probability 
for onc-photon input by pi{t). We plot pi{t) in Fig. [US')! where the signal frequency is set at uJsI^tt = 10.007 GHz 
[the upper impedance-matching spot in Fig. Hd)]. It is observed that pi{t) increases rapidly during the signal pulse 
duration (|t| < 1/2): the excitation probability roughly evolves as pi{t) ~ / (ir|/s(T)p. At the end of the capture 
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time (ns) time (ns) 


FIG. 5: Time evolution of the qubit excitation probability, for a signal pulse containing zero photon: the actual evolution po(f) 
(red solid) and the adiabatic approximation Po(^) (blue dashed). The parameters used are: Q.dl'in = fl™^/27r = 13.2 MHz, 
I = 100 ns, /3 = 2, and tu = 20 ns in (a) and tn = 30 ns in (b). 



time (ns) pulse length (ns) (MHz) 

FIG. 6: (a) Time evolution of the qubit excitation probability pi{t) for a signal pulse containing one photon (red solid). The 

result for zero-photon case, Po{t), is also plotted for reference (blue dashed). The parameters used are: I = 100 ns, /3 = 2, 
w = 30 ns, and iJsl2n = 10.007 GHz. (b) Dependence of the detection probability Pi on the signal pulse length I for various 
values of 7 / 271 : 0.1 MHz (red solid), 0.02 MHz (blue dashed) and 0 MHz (green dotted), (c) Dependence of the detection 
probability Pi on the drive amplitude fid and the signal photon frequency lJs- 


stage, a high detection probability Pi = pi{tf) close to unity is achieved, in clear contrast with the zero-photon case 
where Pq is negligibly small. 

Figure [HKb) shows the dependence of Pi on the signal pulse length 1. Here we assumed three values of the qubit 
decay rate 7 . In the ideal limit of infinite qubit lifetime (7 = 0 ), the detection probability increases monotonically in 
I and a high probability close to unity is attained ior I :s> k~^. This is because a longer pulse is advantageous for the 
destructive interference between the input and elastically scattered photons, which is the key physical mechanism for 
the deterministic switching of the A system. However, in reality, the detection probability is decreased by the decay 
of the qubit into other channels during the pulse duration. A longer pulse is disadvantageous in this regard. We 
thus have an optimum pulse length /opt that maximize the detection probability. Under a practical qubit decay rate 
( 7/277 =0.1 MHz), Pi amounts to 0.89 at /opt ~ 100 ns. If the qubit lifetime is improved five times ( 7/277 =0.02 MHz), 
Pi reaches 0.96 at /opt ^ 200 ns. Note that these detection probabilities represent the switching probability in the 
capture stage of Fig. [3] which does not include the potential errors in the qubit readout stage. 

Figure [He) plots Pi as a function of the drive power and the signal frequency for a signal pulse length / = 100 ns. 
We observe qualitative agreement with the reflectivity plot of Fig. |Hd): Pi is maximized when the drive amplitude 
satisfies the impedance-matching condition (H^ = H™^) and the photon frequency is tuned closely to wai or W 41 . 


C. Two-photon input 

As we discussed in previous subsections, our device discriminates with a high efficiency whether the signal pulse 
contains a photon or not. Here we present the results for two-photon input to observe the multiphoton effects in 
the signal pulse. In Fig. IZKa), we compare the time evolution of pi{t) and P 2 {t)- When the signal photons are out 
of resonance with the A system ( 07^/277 = 10.030 GHz) and therefore the excitation probabilities are low, P 2 {t) is 
roughly twice as large as pi{t). The two photons excites the qubit independently in this case. In contrast, when the 
signal photons are on resonance ( 07^/277 = 10.007 GHz) and a high efficiency is attained for the one-photon input. 
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FIG. 7: (a) Time evolution of the qubit excitation probabilities, pi{t) (red solid) and P 2 {t) (blue dashed). The input photons 

are on resonance with the A system for — 10.007 GHz and are out of resonance for — 10.030 GHz. (b) Ratio of 

the detection probabilities, P 2 /-P 1 , as a function of fid and ujs. 
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FIG. 8: Dependence of the detection probability on the pulse length I for various pulse shapes: Gaussian (red solid), square 
(blue dashed), and exponential (green dotted). The qubit decay rate 7/27r is 0 MHz in (a), 0.02 MHz in (b), and 0.1 MHz in 
(c). 


P 2 {t) becomes comparable as pi{t) due to the saturation of the A system. Figure [T^b) plots the ratio of the detection 
probabilities, P 2 /P 1 , as a function of and uJs- By comparing Figs.lHlJc) and[71[b), we observe that P 2 — 2Pi holds 
when Pi < 0.2 and P 2 < Pi holds when Pi > 0.7. 


D. Pulse-shape dependence 


Up to here, we assumed a Gaussian pulse for the signal photon. Here, in order to observe the effects of the pulse 
shape, we consider square and exponential shapes for the signal pulse: 


/rw = 


Vi ' 

frnt) = 

/21n2 

V / ^ 


— iiOgt 


(|t| < 1 / 2 ) 

m > V2) ’ 

r 2-{t/i+i3l2) 

X < 

10 


{t > -pi/2) 
it < -/tl/2) ■ 


(17) 

(18) 


Note that the amplitude decays exponentially in f™'^, as in the single photons generated by spontaneous emission of 
an emitter. Regarding the drive pulse of Eq. m, we set /3 = 1 (/3 = 3) for the square (exponential) shape to cover 
the signal pulse. We employ the same definition for the detection probability: Pi = pi{tf), where tf = /31/2 + 50 ns. 

Figure [8] shows the dependence of the detection probability on the pulse length 1. For the cases of infinite qubit 
lifetime [Fig. |8l[a)], the detection efficiencies increase monotonically. A high detection probability exceeding 0.9 is 
obtained for I > 100 ns, regardless of the pulse shape. However, if we take account of the finite qubit lifetime [Figs.jSJb) 
and (c)], the detection probability decreases due to the qubit decay. A square pulse is advantageous in this regard, 
since we can use a shorter drive pulse (smaller /3) and suppress the qubit relaxation during the capture stage. Overall, 
regardless of the pulse shape, a high detection efficiency exceeding 0.9 is possible if the qubit has a sufficient lifetime 
(7/27r < 0.02 MHz). 
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FIG. 9: (a) Time evolution of the qubit excitation probability Pg{t) in the reset stage, starting from |(;,0). The mean photon 

number in the reset pulse (n) is 0 (red solid), 10 (blue dashed) and 20 (green dotted). The three lines are mostly overlapping. 
The parameters used are: Qdl2n — 44 MHz, = 9.860 GHz, I = 100 ns, /3 = 2, and w = 30 ns. (b) Time evolution of 

Pe{t), starting from |e, 0). The natural qubit decay is also plotted by the thin line, (c) Qubit excitation probability for (n) = 10 
at the end of reset stage starting from |e, 0). 


VI. RESET OF THE SYSTEM 

In the present device, the qubit-resonator system is automatically reset to the ground state \g, 0) through the qubit 
decay. However, in order to shorten the dead time of the detector, we reset the system through microwave transitions. 
In the reset stage of Fig. |31 we simultaneously apply a drive pulse to the qubit through WG2 and a reset pulse to the 
resonator through WGl: the drive pulse sets the system to the A mode, and the reset pulse induces the |2) — \u) —>■ |1) 
transition (rt = 3 or 4). The drive pulse profile is the same as that in the capture stage, fd{t) of Eq. (03 . However, 
as we observe below, a stronger drive (Q^ > is advantageous in this stage. As the reset pulse, we use a classical 

Gaussian pulse. 


Mt) = ^ X ^ (19) 

where (n) is the mean photon number in the reset pulse. Hereafter, we set uJr close to W32 and use |2) —>■ |3) —>■ |1) 
transition for resetting. Furthermore, we fix the pulse length at I = 100 ns and /3 = 2, and set the initial and final 
moments of the reset stage at ti = —131/2 — 50 ns and tf = [31/2 + 50 ns. 

In Fig. IHKa), we show the time evolution of the qubit excitation probability Pg{t) in the reset stage, starting from 
|g, 0). For ti < t < —131/2, an adiabatic |g, 0) —>■ |I) transition is induced by the drive pulse. For —131/2 <t < 131/2, 
the excitation probability is perturbed by the reset pulse onl;^ slightly, which is due to the large detuning of lOr from 
both W31 and W41. For (31/2 < t < tf, through an adiabatic |1) — >■ |(;,0) transition the system returns to the ground 
state as the drive pulse is switched off. At the end of the stage, the unwanted qubit excitation is suppressed nearly 
completely. For example, the excitation probability Pg =Pg{tf) is below 0.005 for (n) < 20. 

In Fig. [91(b), we show the qubit excitation probability Peit) starting from |e, 0). For ti < t < —/31/2, an adiabatic 
|e, 0) —|2) transition is induced by the drive pulse. For —(31/2 < t < 131/2, we observe rapid decrease of Pe(i), 
which is due to the |2) —> |3) —>■ |I) transition induced by the reset pulse. The reset pulse power (n) determines the 
transition rate. When (n) > 10, the system transits to |I) nearly completely at < = f3l/2. For (31/2 <t <tj, through 
an adiabatic |I) —?► |5,0) transition the system is reset to the ground state. The excitation probability Pe at the end 
of the stage is 0.015 (0.007) for (n) = 10 (20). This is in clear contrast with the slow natural decay of the qubit, 
where the excitation probability remains 0.83 at t = tf. 

In Fig. EKc), Pe is plotted as a function of the drive amplitude Qd and the reset pulse frequency We observe 
that the impedance-matching condition (fl^ = is not necessarily required in the reset stage. The system is reset 

efficiently when cUr — W32 and which implies K31 > 1^32 from Fig. |2](a). Although the reset efficiency per 

one reset photon is lowered by the unbalanced decay rate, many photons ((n) > 10) involved in the pulse enables 
highly efficient resetting. 
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VII. SUMMARY 

We theoretically analyzed the microwave response of a driven qubit-resonator system and demonstrated its excellent 
performance as a detector of itinerant microwave photons. The detection of microwave photon proceeds by switching 
two modes adiabatically: Under a proper qubit drive, the system functions as an impedance-matched A system, which 
efficiently absorbs itinerant photons (A mode). Due to the destructive interference between the incident and elastically 
scattered photons, a high efficiency close to unity is attained for temporally long pulses, regardless of their shapes. 
In contrast, without a qubit drive, the system functions as an oscillator whose resonance frequency depends on the 
qubit state (I mode). We can perform dispersive readout of the qubit state in this mode. The detector is operated 
in the following cycle: In the capture stage, we set the system in the A mode and map the presence (absence) of 
a photon in the signal pulse to the excited (ground) state of the qubit. In the readout stage, we set the system in 
the I mode and measure the qubit state. In the reset stage, in order to shorten the dead time of the detector, we 
restore the system to its ground state through microwave transitions. The detection efficiency readily exceeds 0.9 for 
realistic parameters and the dark counts are suppressed nearly completely. The present scheme provides a long-sought 
single-photon detector in the microwave domain and widens the options in microwave quantum-optics experiments. 
It could be used for the characterization of photon number statistics in nonclassical itinerant microwave fields and 
the measurement of flying qubits encoded in propagating microwave photons. 
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Appendix A: Relation between coherent-state and Fock-state inputs 


Here we derive Eq. m, which connects the expectation values for a coherent-state input to those for a Fock-state 
input [S^. We denote a coherent state by |a) and a Fock state by \m) (m = 0,1, • • •). We denote the expectation 
value of an operator S{t) for the coherent-state input by {S)c = (a|S'|a) and its component proportional to (a*)™a" 
by (5')™"'. On the other hand, we denote the matrix elements in the Fock bases by {S)mn = (m|S'|n). 

Using the Fock-state expansion of a coherent state, |q;) = (‘^)c written as 


(5)c = e 






{S)r 


(Al) 


Picking up the terms proportional to (Q;*)°a°, (a*)^a^ and (a*)^a^, we have (S')™ = (S)oo, {S)l^ = (S)ii — (S)oo 
and {S)f = {8)22/2 - {S)n + {S)oo/2. Therefore, 

(S)oo = (S)™, (A2) 

(S)n = {S)/° + {S)l\ (A3) 

(S)22 = (S)™ + 2(S)“ + 2(S)f. (A4) 
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